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Abstract 

The centers of elliptical galaxies host supermassive black holes that signif- 
icantly affect the surrounding interstellar medium through feedback resulting 
from the accretion process. The evolution of this gas and of the nuclear emis- 
sion during the galaxies' lifetime has been studied recently with high-resolution 
hydrodynamical simulations. These included gas cooling and heating specific 
for an average AGN spectral energy distribution, a radiative efficiency declining 
at low mass accretion rates, and mechanical coupling between the hot gas and 
AGN winds. Here we present a short summary of the observational properties 
resulting from the simulations, focussing on 1) the nuclear luminosity; 2) the 
global luminosity and temperature of the hot gas; 3) its temperature profile and 
X-ray brightness profile. These properties are compared with those of galaxies 
of the local universe, pointing out the successes of the adopted feedback and the 
needs for new input in the simulations. 
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1. Introduction 

Supermassive black holes (MBHs) at the centers of bulges and elliptical 
galaxies played an important role during galaxy formation and evolution, as 
revealed by remarkable correlations between their masses and some host galaxy 
properties (e.g., Ferrarese & Merritt 2000, Gebhardt et al. 2000) and by many 
studies (e.g., Merloni et al. 2004; Sazonov et al. 2005; Di Matteo, Springel & 
Hernquist 2005, Hopkins et al. 2006). An important aspect of the coevolution 
process is the interaction between the energy emitted by the accreting MBH and 
a galactic interstellar medium (ISM) continuously replenished by normal stellar 
mass losses. In the absence of feedback from a central MBH, this ISM develops 
a flow directed towards the galactic center, accreting >1Mq yr~^ in a process 
similar to a "cooling flow" (e.g., Ciotti et al. 1991, David et al. 1991, Pellegrini 
& Ciotti 1998). Therefore a few basic questions arise: 1) what is the fate of the 
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large amounts of gas accreting towards the center during the galaxies' lifetime, 
and not observed at any wavelength (e.g., Peterson & Fabian 2006)? Just ~1% 
of the mass made available by stars is in the mass of present epoch MBHs 
(Ciotti & Ostriker 2007). 2) how much radiative and mechanical energy output 
from the MBH can effectively interact with the surrounding ISM? 3) what are 
the resultant masses of the MBHs at the present epoch? and finally, 4) why 
bright AGNs arc; not common in the spheroids of the local Universe, given the 
expected mass accretion rate (e.g., Fabian & Canizares 1988, Di Matteo et al. 
2003, Pellegrini 2005a)? 

Recently, the interaction of the MBH with the inflowing gas has been stud- 
ied with high resolution ID hydrodynamical simulations by Ciotti & Ostriker 
(2007), Ciotti, Ostriker & Proga (2009, hereafter paper I) and Ciotti & Ostriker 
(2009, hereafter paper H). These simulations included a detailed and physically 
based treatment of the radiative energy output from the MBH and its trans- 
fer in the ISM, and of the mechanical energy from AGN winds. They showed 
recurrent brief flaring of the nucleus, with the galaxy seen alternately as an 
AGN/star burst for a small fraction of the time and as a "normal" quiescent 
elliptical for much longer intervals. Accretion feedback proved effective in sup- 
pressing long lasting cooling flows and in maintaining MBH masses within the 
range observed today (since gas is mostly lost in outflows or starbursts); a ma- 
jor role in regulating the flow evolution was played also by type la supernovae 
(hereafter SNIa's). 

Here we present a preliminary investigation of the spectral and morpholog- 
ical appearance of this class of models in the X-ray band. In particular, we 
focus on the hot ISM and the nucleus during the galaxy's lifetime, both in 
quiescence and during outbursts of activity, and we compare their properties 
with the observations collected for galaxies of the local universe by ROSAT 
and more recently by Chandra. Section 2 briefly describes the chosen repre- 
sentative model, Section 3 its observational properties and a comparison with 
observational results. Section 4 summarizes the conclusions. We anticipate here 
some general results for this class of models, and a more extensive description 
will be given in Pellegrini, Ciotti & Ostriker (2009, hereafter paper III). 

2. The representative model 

The main properties of feedback modulated galactic accretion flows are given 
in Ciotti & Ostriker (2007) and recently in papers I and II, together with a full 
description of the code used for the simulations. Briefly, the code integrates 
the time-dependent ID Eulerian equations of hydrodynamics, with a logarith- 
mically spaced and staggered radial grid, extending from 2.5pc from the central 
MBH to 250 kpc. The code calculates self-consistenly the source and sink terms 
of mass, momentum and energy, associated with the evolving stellar population 
(stellar mass losses and SNIa's events), nuclear starbursts, accretion and MBH 
feedback. Gas heating and cooling are calculated for a photoionized plasma 
in equilibrium with an average quasar spectral energy distribution (Sazonov et 
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al. 2005), the resulting radiation pressure and absorption/emission are com- 
puted and distributed over the ISM from numerical integration of the radiative 
transport equation. 

Of interest here is that the adopted radiative efficiency of material accreting 
on the MBH at a rate M is e = 0.1 x 100m/(l + 100m), where m = M/MEdd 
is the Eddington-scaled accretion rate, so that e ~ 0.1 at large mass accre- 
tion rates (when m ^ 0.01), and declines in a RIAF-like fashion e ^ lOm for 
m<0.01 (as for radiatively inefficient accretion flows, Narayan & Yi 1994). The 
mechanical feedback implemented is that of quasar outflows [e.g., Chartas et al. 
(2003), Crenshaw ct al. (2003)], as modelled numerically by Proga (2003). The 
mechanical efficiency scales with the Eddington ratio {I = L/Lsdd), reaching a 
maximum value for Z = 2 of 3 x 10~^ to 5 x 10~^ in different models. Note that 
only the mechanical feedback due to the Broad Line Region wind is activated 
(the nuclear jet properties will be inserted in a future work). 

Here we consider a representative model from the latest set of simulations 
(paper II) , for an isolated typical elliptical galaxy with a central stellar velocity 
dispersion a = 260 km s~^, a B-band luminosity Lb = 5 x IO^^Lb,©, a stellar 
mass described by the Jaffe law with a mass-to- light ratio M^/Lb = 5.8, an 
effective radius = 6.9 kpc, a dark matter halo with equal amount of dark 
and visible mass within i?g, a standard SNIa's rate; the maximum value of the 
mechanical efficiency is set to 10~^. The simulations begin at a galaxy age of ^ 2 
Gyr (i.e., a rcdshift z ^ 2, the exact value depending on the epoch of elliptical 
galaxy formation, usually put at z>2) and soon AGN outbursts develop, each 
followed by an abrupt drop of the accretion rate, and by times during which the 
galaxy is replenished of gas by the stellar mass losses. 

The typical behavior of the gas during an outburst starts with the off-center 
growth of a shell of denser gas (at a radius of ~ 0.5 — 1 kpc) that progressively 
cools, falls and reaches the center (i.e., a cooling collapse occurs). Then a 
radiative shock quickly [in ~ (1 — 2) x 10^ yrs] produces an outward moving 
shell of cold and dense gas, that slows down until it falls back towards the center, 
accumulating far more material than the first shell, and giving origin to a larger 
accretion and feedback episode. After some more events of this kind, a major 
shock leaves behind a very hot and dense center, and causes substantial galaxy 
degassing. The gas then cools and resumes its subsonic velocity, the density 
starts increasing again and the cycle repeats. Overall, during outbursts both 
hotter and colder regions are continuously created, the hotter ones mainly at 
the center (due to gas compressed by the falling shells, or to shocks) and the 
colder ones due to the gas in the falling shells or in the radiative shocks. 

We remark that the representative model has fairly standard input parame- 
ter values, and its general properties described below are typical of this class of 
models. A large set of runs exploring the impact of different values for the feed- 
back parameters (e.g., the wind opening angle, the peak mechanical efficiency 
and its dependence on I) showed present-day properties similar to those of the 
adopted representative model (papers I and II). 
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3. Comparison of simulated and observed X-ray properties 

We consider in turn the evolution of the luminosity of the nucleus (Sec- 
tion 13. ip . of the total luminosity and emission weighted temperature of the 
hot gas (Section 13. 2p . of the temperature profile (Section 13. Sp and of the sur- 
face brightness profile fSection 13. 4p . The emission is calculated over the full 
Chandra sensitivity band (0.3-8 keV), and in two separate bands, 0.3-2 keV 
and 2-8 keV, by means of the APEC code within the package XSPEC for the 
X-ray data analysis (see paper III for more details). 

Since in the simulations the treatment of feedback is physically based, not ar- 
bitrarily chosen and tuned to reproduce observations, any agreement or discrep- 
ancy of the resulting model properties with observations is relevant to improve 
our understanding of the MBH-ISM coevolution. 

3.1. The nuclear luminosity 

The time evolution of the bolometric nuclear luminosity shows strong in- 
termittencies at an earlier epoch, reaching the Eddington value, and becoming 
rarer and rarer with time, until a smooth, very sub-Eddington phase establishes 
(Fig.[T]). This behavior follows from the secular decrease of the stellar mass loss 
rate (Fig. [T|) , which produces longer and longer times for the replenishment of 
the galaxy with gas. Towards the present epoch (a galaxy age of ^ 12 Gyr) 
the mass accretion rate is M O.OlMo/yr, that is m ~ 1.1 x lO""^, therefore 
accretion has entered the RIAF regime and e ~ 0.01 (Sect. [2]). The nuclear 
bolometric luminosity is Lboi = 2 x lO'*'^ erg s^^ and its Z = 2 x 10^**, in agree- 
ment with the fact that in the local universe the fraction of MBHs approaching 
their Eddington limit is negligible. For example, in the statistically complete 
Palomar spectroscopic survey of 486 nearby galaxies only ~ 50% of ellipticals 
show emission line nuclei, mostly of low level {Lhu < 10^" erg/s); for LINERs, 
a sub-sample largely dominated by ellipticals, Lboi ranges from 10^^ to 10"*^ erg 
s""'^, with a peak at L^oi ^ 10^^ erg s~^, and the mediarQ I = (0.5 — 1) x 10^^ 
(Ho 2008). Seyferts excluded, nuclei of all types have Li,oi < 10^"^ erg s^^. 
Therefore the representative model may need a reduction of its L^oi , in addition 
to that produced by the adopted low radiative efficiency of RIAFs, to be more 
representative of the typical low luminosity nucleus of the local universe. 

Indeed, many nuclei of nearby ellipticals studied in detail in the X-ray band 
show emission at an extremely low level. Their 0.3-10 keV luminosities Lx,nuc, 
for a sample within ~ 50 Mpc, with or without radio emission and residing in 
all kinds of environments, range from fewxlO'^* erg s~^ to 10^^ erg s~^, with 
Lx.nuc/ LEdd as low as 10^^ — 10^® the most common (Pellegrini 2005a, b; see 
also Soria et al. 2006; Gallo et al. 2008). Lx,nuc of the model at the present 
epoch can be derived from Lf,oi adopting a correction factor appropriate for the 
spectral energy distribution of a RIAF, that is <0.2 for low luminosity MBHs 



^Ljjo! derives from the nuclear 2-10 keV emission, with Lhai/L^—io keV = 16 as appropriate 
for low luminosity AGNs, and L^dd from MBH masses estimated from the MBH-cr relation. 
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(Mahadevan 1997). This gives Lx,™ic^4 x 10''^ erg s~-^, and Lx,nuc/ Lsdci^^ x 
10~^, both values again quite higher than typically observed. 

A reduction of the model nuclear luminosity could be accomplished by in- 
serting the mechanical feedbaclH of a jet or of a nuclear wind, that should 
naturally develop in low radiative efficiency accretion (Blandford & Begelman 
1999), and should reduce further the mass accretion rate (e.g., Di Matteo et al. 
2003). In fact, the Bondi accretion rate for the hot ISM around many nearby 
MBHs, coupled to the RIAF radiative efficiency, corresponds to luminosities 
still higher than observed, indicating the need for a reduction of the mass that 
actually reaches the MBH (Loewenstein et al. 2001, Pellegrini 2005a, Pellegrini 
et al. 2007a, b). Another indication of the relevance of mechanical feedback is 
the observed link between the radio luminosity and the X-ray gas asymmetry 
index for ~ 50 nearby ellipticals (Diehl & Statler 2008a). 

Finally, note that Lx,nux is observed to be independent of the MBH mass and 
of the Bondi accretion rate (Pellegrini 2005a). Likely, neither Bondi accretion 
nor simple RIAFs realistically describe the interaction between the outgoing 
energy flux and the incoming mass flux at the nucleus, while the fluctuations 
typical of feedback cycles can account for the observed lack of relationship. 

3.2. Global gas luminosity and temperature 

The time evolution of the gas emission shows peaks during outbursts, that 
become broader and less high with time increasing (Fig. Soft and hard 
emission oscillate in phase and present the same overall behavior, with the hard 
one keeping <100 times lower. Hard emission during outbursts would be difficult 
to distinguish from the contribution of unresolved binaries even with Chandra, 
if extended (e.g., Trinchieri et al. 2008), but it could be detected if centrally 
concentrated (see also Sect. 13. 4| ). 

The largest compilation of the global X-ray emission of early type galaxies 
of the local universe is based on ROSAT observations of a relatively unbiased 
sample of 401 objects, homogeneously analyzed (O'SuUivan et al. 2001). The 
best fit correlation with the B-band luminosity of the emission mostly due to hot 
gas (Lxjsm), after subtraction of the stellar contribution, predicts Lxjsm ~ 
10^"'^ erg s""'^ for Lb = 5 x 10^^ Lb,q as the model galaxy. However at this 
Lb, due to the large scatter around the best fit, the observed Lx.ism ranges 
from 10'*'^ erg s~^ up to few xlO^^ erg s~^ (though Lx.ism > fewxlO^^ erg 
belong to galaxies with an important contribution from the intragroup or 
intracluster medium). At an age of ^ 10— 12 Gyr the model has Lxjsm ^ 10*^° 
erg (Fig. 12), on the lower end of the observed range. This indicates that 
degassing is too efficient in the simulations, at least at the present epoch. This 
problem would be alleviated by an external pressure from an outer medium 
(e.g., Vedder et al. 1988), a more appropriate "boundary" in fact for ellipticals 
that usually are not isolated (it will be considered in future simulations). 



■^We recall that the simulations include only the mechanical feedback due to the Broad 
Line Region wind (see Sect. [2l and Fig. 1 in paper I). 
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The lower panel of Fig. [5] shows the evolution of the soft (< T >s) and 
hard (< T >h) emission weighted temperatures within the optical i?e, and their 
complex behavior during outbursts. The sharp and high peaks in < T >;i 
correspond to the onset of very hot regions, while the decrements in < T >s 
are due to a dense cold shell preceeding the major burst and to cold gas swept 
by the radiative shocks produced by the outburst (see Sect. [5]). The 0.3-8 keV 
emission weighted temperature (not shown in Fig. [2]) is almost coincident with 
< T >s, except for those short times with a very hot gas component. 

ROSAT observations showed a correlation, though with a large degree of 
scatter, between the temperature weighted with the soft emission of the whole 
galaxy (T) and the central cr (O'SuUivan et al. 2003), predicting T ~ 0.7 keV 
(with observed values from ^ 0.4 to ^ 1 kcV) for a — 260 km s^^ as for the 
model galaxy. At the present epoch the model has T ~ 0.5 keV, lower than the 
best fit prediction but within the observed range, and < T >s= 0.7 keV (Fig.[2|). 
Likely the observed T refers to a region less extended radially (and typically also 
hotter) than used for the model galaxy. Average temperatures within Re have 
been obtained with Chandra (Athey 2007), whose high angular resolution allows 
for a better subtraction of the point (stellar) source contribution and then gas 
temperatures closer to the true ones. For 22 Athey's galaxies with log Lb = 10-5 
to 10.8, similar to that of the model, the average 0.3-8 keV emission weighted 
temperature is 0.60 keV, close to < T >s though shghtly lower. 

Summarizing, the gas content of the model is lower than the average observed 
for a galaxy of its size; this shows the need for a lower SNIa heating, or a 
higher/more concentrated gravitating mass, or the confining effect of an external 
medium. The average temperatures instead agree reasonably with observations, 
and modifications to the model should not increase the average temperature 
within i?e, while they could increase a little that over the whole galaxy. 

3.3. Temperature profiles 

The emission-weighted and projected temperature profiles are shown during 
quiescent inter-burst periods in Fig. [31 and during the last outburst at 7.5 
Gyr in Fig. 21 During quiescent times the profiles are smooth, with the tem- 
perature monothonically decreasing for increasing radius (from 1 keV at a 
radius of ~ 100 pc to ~ 0.4 — 0.5 keV at ~ 20 kpc). Negative temperature gra- 
dients are typical of inflowing gas in steep potentials (e.g., Pellegrini & Ciotti 
1998), just due to compressional heating. Diehl & Statler (2008b) suggest that 
negative temperature gradients could be the sign of localized heating by a cen- 
tral, weak AGN. In fact in the models, even during quiescent times, accretion 
is present (though with a low radiative efficiency, Sects. [H and 13. 1|) , and this 
has an additional heating effect over gravitational compression on the gas at 
the center. Negative gradients are commonly observed among ellipticals, as re- 
vealed by Chandra for an increasing number of galaxies (Fukazawa et al. 2006, 
Athey 2007, Diehl & Statler 2008b), and the temperature roughly halves from 
the center to the galaxy outskirts, as in the model. A few galaxies have central 
temperatures of ~ 1 kcV (as in Fig. [3]) , while most show temperatures between 
0.5 and 1 keV at their innermost radial bin. However, the model temperature 
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calculated as an average for a central bin extending out to 0.5-a few kpc, as for 
observed galaxies, will be lower than 1 keV. While a more detailed comparison 
is postponed to paper III, we conclude here that no additional heating seems to 
be required in the central region (as found in Sect. 13. 2p . 

In other observed profiles the temperature increases outward, or keeps roughly 
flat, or is "hybrid", firstly decreasing from the central value to a minimum and 
then increasing out to the edge of the galaxy (Diehl & Statler 2008b). These 
profiles could be produced in the model by environmental effects currently not 
included (e.g., Vedder et al. 1988, Mathews & Brighenti 2003). A temperature 
decreasing towards the center is also observed, but is never shown by the model 
during quiescence. A central drop in temperature and a shape resembling the 
hybrid profile are shown during particular phases of outbursts, when the major 
shock is moving outwards or the cold shell is forming (see Fig. |4] that illus- 
trates the messy temperature distribution during the 0.2 Gyr of an outburst, 
consequent to the gas flow behavior described in Sect.[2|). 

We noted (Sect. 13. ip that the stationary hot accretion phase needs a reduc- 
tion of the accretion rate, as can be provided by a jet. In a recent Chandra 
analysis the temperature profiles outside ^ Re showed gradients that switch 
monothonically from positive to negative going from high to low radio luminosi- 
ties within three Rf. (Diehl & Statler 2008b). To be consistent with this result, 
the jet addition to the simulations should heat the gas outside ^ Re ■ 

Finally note that the temperature profile of the elliptical NGC3411 shows 
an intriguing off-center dip (O'SuUivan et al. 2007) that is highly unusual for 
standard cooling flow models but is present in a few profiles in Fig. ID This dip 
corresponds to a ripple in the observed X-ray brightness profile, that is also a 
frequent feature in the model (Sect. [3^ . 

3.4- Brightness profiles 

The evolution of the surface brightness profile is shown in Fig.[5l right before 
and during the last major outburst, and in Fig. [51 during the subsequent slow 
evolution until the present epoch quiescent state. In Fig.[5]curves labelled 1 and 
2 show two subsequent times of the off-center shell formation, lasting ~ 0.5 Gyr. 
Curve labelled 3 shows a shock moving outwards after the major outburst (in the 
upper panel) and the presence of very hot gas at the center (lower panel); this 
phase lasts for 2 x lO'' yr. Curve labelled 4 shows the result of the degassing 
caused by the passage of the shock waves: the galaxy has a low gas density 
and subsonic perturbations remain at a radius of a few tens of kpc; the time is 
<0.2 Gyr since the outburst started. Note how the main features in the profile 
are above the levels of background and of unresolved stellar emission (described 
in the caption of the figures). The high-temperature, high-density phase at 
the center (curve 3) is very brief, and unlikely to be observed. Disturbances as 
shells and ripples lasting <0.2 Gyr are more likely to be observed. In fact, many 
nearby galaxies show these characteristics in their images, as recently proven 
with Chandra studies (Diehl & Statler 2008a, Machacek et al. 2006). 

Interestingly, three bright ellipticals imaged with Chandra (Loewenstein et 
al. 2001) all show a flatfish X-ray brightness profile within the central 3-10 
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arcsec (<1 kpc), that is not possible to reproduce with inflow models, while it 
resembles the profile of the "pre-burst" phase (curves 1 and 2 in Fig. [5]). 

4. Summary and conclusions 

We have calculated the X-ray properties of a galaxy model representative of 
a recent investigation of feedback modulated accretion including the combined 
heating of radiation and AGN winds (papers I and II), and we have compared 
them with observations of the local universe. Besides important achievements 
as limiting the growth of the MBH mass to observed values (papers I and II), 
thanks to the present analysis these models also reveal important aspects of the 
MBH-ISM coevolution, as detailed below. 

At the present epoch a highly sub-Eddington {I ^ 10^'*) phase establishes 
with accretion in the low radiative efficiency regime. Though within the range 
observed, the nuclear emission is quite larger than the most frequently observed 
values [/ (0.5 — 1)10^''']. This suggests that an additional mechanism reducing 
further the mass available for accretion is important, as could be provided by a 
nuclear jet or wind from a RIAF. 

During the many outbursts, the gas emission Lx.ism presents peaks that 
become broader as time increases; sharp peaks and decrements are shown re- 
spectively by the hard and soft emission weighted temperature. At the final 
quiescent times, Lxjsm is on the lower side of those observed, which suggests 
degassing may be less efficient, or the confining agents (gravitating mass, ex- 
ternal medium) more efficient. The emission weighted temperatures instead 
compare well with observations. 

The profiles of brightness and temperature (decreasing outward) in quies- 
cence resemble those of many local galaxies. Outbursts produce features in the 
brightness profile that are detectable with Chandra, last for <0.2 Gyr, and 
could match part of the widespread gas disturbances observed in local galaxies. 
The hot bubbles inflated at the center produce central peaks in the hard band 
brightness profile, but their short duration makes them difficult to observe. 

In conclusion, a form of feedback lowering further the nuclear luminosity 
during the stationary hot accretion phase seems important. This should not 
increase the central temperature that is already consistent with observed values. 
A jet interacting mainly outside ^ 0.5 — IRe would then be preferable; it should 
also heat the gas there, to bring the overall temperature profile in accord with 
those observed (increasing outward) at high radio luminosities. The additional 
feedback should not lower Lxjsm- 

Finally, we comment on how Lxjsm and the gas temperature might scale 
with galaxy properties, and with how much scatter. A change of the mass 
content and distribution, SNIa rate and external pressure will have an impact 
on the gas content, and also on the possibility of having outbursts until the 
present epoch, or no outbursts at all. For example, the SNIa heating may 
keep small galaxies in an outflow so that outbursts never happen; due to a 
larger dark matter content, galaxies as the representative model may never 



8 



experience a huge degassing, and may host nuclear outbursts until the present 
epoch. In general, we expect a trend of Lxjsm with Lb mainly determined 
by the galaxy structure and SNIa's heating, as described by older studies (e.g., 
Pellegrini & Ciotti 1998), to which scatter is added by feedback, but only for 
luminous galaxies (as the representative model). The frequency of outbursts at 
the present epoch is difficult to estimate, given the many parameters involved, 
but it should be very small (papers I and II). The hot gas disturbances shown 
by Chandra have been mostly attributed to an ongoing or recent jet activity 
(e.g., Diehl & Statler 2008a,b), though evidence for this is not always present 
(e.g., NGC4552, Machacek et al. 2006). A statistically complete sample of 
galaxies with information on the nature of activity and of the gas disturbances 
will hopefully be built soon from the large Chandra database. 
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Figure 1: Left panel: the evolution of the bolometric luminosity produced by nuclear accre- 
tion. The upper line is the Eddington luminosity (note the slow and overall modest increase 
of the MBH mass). Right panel: the evolution of the mass budget: the stellar mass loss rate 
(green line), the gas mass leaving the galaxy at a radius of 10/?e (black line) and the mass 
accretion rate on the MBH (red line). 
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Figure 2: Time evolution of the total gas luminosity and emission weighted temperature within 
the optical effective radius Re, in the 0.3-2 keV band (red lines) and in the 2—8 keV band 
(blue lines). 
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Figure 3: Projected temperature profiles weighted with the emission in the 0.3-8 keV band, 
during the interburst times t=6.9 Gyr (solid line), 8.1 Gyr (dotted), 9.1 Gyr (dashed) and 
11.1 Gyr (dot-dashed line). 
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Figure 4: Projected temperature profiles for the gas emission during the last outburst (see 
Sect. [2] for more details). Upper panel: the shell is forming at a radius of ~ 1 kpc and 
approaching the center (solid line); a shock is moving outward after the first burst (dotted 
line, after 54 Myr; dashed line, after 4 Myr), creating a major shell that falls back to the 
center (dot-dashed line, after 8 Myr). Central panel (the time step between each line is 2 
Myr): after the shell has reached the center, a shock is moving outward, creating a central 
very hot bubble (solid line, dotted line), that quickly cools (dashed and dot-dashed line). 
Bottom panel: the shock is moving outwards [solid (-1-4 Myr), dotted (+6 Myr), dashed (-1-16 
Myr) lines]; quiescence is resuming (dot-dashed line, -1-0.1 Gyr). The time elapsed from the 
beginning of the outburst (Fig. [J] upper panel, solid line) is 0.2 Gyr. 
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Figure 5: The brightness profile during the last major outburst; curves referring to subsequent 
times are labeled from 1 to 4 (see Sect. 13.41 for a description of the main features). Count 
rates in the y-axis and background level (straight horizontal line) are for the Chandra ACIS-S 
detector. The red line shows the profile of unresolved X-ray binaries, appropriate for their 
spectral shape and a representative observation of a Virgo galaxy with an exposure of 200 
ksec (e.g., Kim et al. 2006). 
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Figure 6: The evolution of the brightness profile after the last major outburst: the quiescent 
state establishes again and is kept until the end of the simulation. The solid, dotted, dashed 
lines correspond to times of 7.7, 8.7 and 10.0 Gyr. Count rates and unresolved stellar emission 
profile (thin smoother line) have been calculated as for Fig. \5\ 
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